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In this article, we develop quantum mechanics upon the framework of the quantum mechanical 
Hamilton-Jacobi theory. We will show, that the Schrodinger point of view and the Hamilton- 
Jacobi point of view are fully equivalent in their description of physical systems, but differ in their 
descriptive manner. As a main result of this, a wave function in Hamilton-Jacobi theory can be 
decomposed into travelling waves in any point in space, not only asymptotically. The well known 
WKB-theory will be show up as a special case of the more general theory, we develop below. By 
the example of the linear potential and the harmonic oscillator we will discuss quantum mechanics 
from the Hamilton-Jacobi point of view. Soft boundary value problems as the connection problem 
can be solved exactly. Quantizised energies and Maslov-indices can be calculated directely without 
orthonormalizing wave- functions. Also, we will focus on trajectory-themes, which, in contrast to 
the Schrodinger point of view, follow naturally from the quantum mechanical action function. 
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It must have been a subconscious path, that guided Schrodinger to the derivation of his famous equation in the 
early days of 1926 [J. Schrodinger started with the classical Hamilton-Jacobi equation, transformed the classical 
action function s by s = — iftln^], performed an adventurous variational calculation on the Hamilton-integral and 
arrived at his celebrated result 



$"(x) + p 2 {x)^(x) = 0, p(x) = ^/2m(E~U{x))/h . (1) 

04 , Several attempts have since been made to relate the quantum mechanics perfectely described by |T]) to quantities, 
that are known from classical mechanics, action and momentum; these attempts are twofold. The first one was taken 
by Madelung @ already in 1926. Madelung introduced the polar wave-function, \&(x) = A(x)exp[iSu(x)/H], where 
Sm(x) is the Madelung action function and both A, Sm(%) are real. Madelungs way was further taken by Bohm 3] 
' in his attempt to build a causal quantum theory. But the initial separation of real and imaginary parts of the wave- 
On, function leads to very complicated equations for amplitude and phase and will not be pursued here. For exaustive 
-i— > 1 information on the Madelung-de-Broglie-Bohm theory see [J, [|[ and for modern application Q . The work of Floyd 0] 
pursued especially the theory of the stationary Hamilton-Jacobi theory with quantum potentials, while Faraggi and 
Matone [8[ probably worked out some fundamental ideas on the geometrical meaning of the quantum potential. In 
the appendix, we will comment especially on '7;] in the light of our results. 

The second way is to initially keep the real and imaginary parts of the quantum action function together, ^>(x) = 
exp[iS(x)/K], leading to the Hamilton-Jacobi equation of quantum mechanics, 

(S'(x)f-ihS"(x)=p 2 (x) . (2) 

Also in 1926, Wentzel, Kramers and Brillouin [9(, starting from ([2]), proposed an expansion of the full quantum action 
function S(x, H) by small values of h, which finally leads to the classical action, pertubed by small corrections of 
order h. The WKB-theory works surprisingly well, but it breaks down in the vincity of quantal regions, where the 
WKB-wave function is plagued by singularities. For an exhaustive review on this topic see [l(| and references therein. 

Various methods, most of them are reviewed in [ll[, have been developed to improve the WKB-theory. The 
most elaborated among these are the Froman-Froman series [l2| and the use of the quasi-linearization method by 
Raghunathan and Vasudevan [l3T ]. The method of quasi-linearization has the advantage, that it initially does not rely 
on the existence of a small expansion parameter, h in our case. Quasi-linearization, as applied to quantum mechanics 
in [l3| is the starting point of our article. We first develop the relations between the Schrodinger theory and the 
Hamilton-Jacobi theory of quantum mechanics, which, uninterpreted, can already be found in [l|. Then, we are going 
to develop the use of quasi-linearization techniques from a systematic and physically reasonable point of view, which 
was neglected in [X3( ] . The systematic use of this technique will lead us to the definition of a quantum correction 
function. This is a function, that describes all quantum effects of a physical system. Further, the properties of the 
quantum correction function will be analyzed in detail, especially their relation to the WKB-series expansion. After 
the properties of the quantum correction functions are clarified, we will apply them to the eigenvalue-problem and 
the question of quantum trajectories. 

In our article, we will develop a quantum theory from the Hamilton-Jacobi point of view. In contrast to the de 
Broglie-Bohm ansatz, we will not initially separate the amplitude and phase of the wave-function. 
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Structural Notes 



In this subsection, we will examine the interrelation between the Schrodinger and Hamilton- Jacobi equation from 
a structural point of view. This will clarify, that the solutions of both equations indeed are related to each other by 
the well known relations 



if?(x) — exp 



, S(x) = -ih\n[9(x)] . (3) 

By introducing a local gradient field, which can easily be identified with the local quantum momentum, 

S(x) = J P{x)dx P(x) = S'(x) , (4) 

the Hamilton-Jacobi equation @ can be cast into 

P 2 {x) -ihP'(x) = p 2 {x) . (5) 
The full solution of the Schrodinger equation (jT]) is a superposition of two linearly independent solutions 

*(ar) = ci*i(x) + c 2 ^ 2 {x) . (6) 

In the language of action, this reads 



^(x) — c\ exp 



c 2 exp 



- h S 2 {x) 



(7) 



We want to emphasize here, that ([7]) allows us to decompose a Schrodinger wave-function into running waves in any 
point in space, not only asymptotically. The statistical properties of running waves are given by 



Pl (x) = 9;(x)9i(x) = exp[-2&Si(x)/h], ji(x) = ±&b i (xy*' i (x) = ±W 3 i {x)p i (x) 



(8) 



Any of the two linearly independent wave functions solve the Schrodinger equation seperately, and any of the two 
action functions solve the Hamilton-Jacobi equation seperately. But when we take the logarithm, we find 



ihln[^{x)} = S(x) 



Ci exp 



S x {x) + S 2 (x) 



-Mx) 



+ c 2 exp 
2 cos 



S 2 (x) 



Si(x) - S 2 (x) 
2h 2 



(9) 



where the Cj's have been absorbed in the phase <f>. The function S(x) is again a solution of the Hamilton-Jacobi 
equation, but Si(x) + S 2 (x) is not. This fact is due to the nonlincarity of the Hamilton-Jacobi equation and its 
physical interpretation is interference. Usually, the two linearly independent solutions ^i(x), Psi 2 (x) will be complex 
conjugates, such that S 2 (x) = —S*(x). The superponed wave function, up to an amplitude factor, then reads 



~ exp [-%{Si(x)}/h] cos 



R{ Sl (x)}/h-^ 



(10) 



II. QUANTUM CORRECTION FUNCTIONS 

To be able to distinguish between the influence of classical mechanics and quantum mechanics on the behaviour of 
a physical system, we propose, that the full quantum momentum field P(x) can be decomposed into its classical part 
and a contribution arising purely by quantum effects. We set 



P(x) = p(x) + Q{x) 



(11) 



The function Q{x) is called a quantum correction function, because it decribes purely quantum effects and thus 
corrects the behaviour of the classical movement, governed by the classical momentum p{x). More precisely, the 
function Q(x) describes the difference between the full quantum momentum P{x) and its classical contribution p(x) 



Q(x)=P(x)-p(x) . 



(12) 
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It can be concluded easily, that Q(x) tends to zero in regions, where the motion is known to be purely classical. 

To derive an analytic expression for the quantum correction function Q(x), we insert ([lip into the Hamiltion-Jacobi 

equation J5]) and find an differential equation for Q{x): 

Q{x){2p{x) + Q(x)} - ih{ P '{x) + Q'(x)} = . (13) 

This nonlinear first order equation can be formally be integrated to yield a selfconsistent formula for the quantum 
correction function Q(x): 

Q(x)= ~exp[~j-f x {2p(y) + Q(y)}dy]J x p'(y)e X p[jrJ y {2p(z) + Q(z)}dz]dy . (14) 

We emphasize, that this expression for Q(x) is exact; so far, no approximation has been made. In the following, we 
will examine the structure of this expression and formulate a systematic perturbational approach. This perturbational 
approach does not, as it is usually done for semiclassical questions, start initially with the expansion in terms of H, it 
will consist of quantum correction functions, that contain every order in h. This will open a more general sight onto 
the behaviour of physical systems near their semiclassical limit. 



A. Perturbation Theory 



The main reason to examine quantum mechanical systems in Hamilton- Jacobi formulation is to clearly distinguish 
between classical and quantum behaviour. As we will see below, the semiclassical limit of Hamilton- Jacobi quantum 
mechanics is open to perturbation theory, while its full quantum description belongs to the non-perturbative limit of 
this theory. The main result of the following subsection is, that in the one-dimensional case, we are able to introduce 
an operator formalism, which is, in its structure, closely related to the wellknown Dirac- or interaction representation. 
The central quantity of this operator formalism is given by 



Q{x) = exp 

where the quantum action function is given by 



-W{x) 



Q(x) exp 



-W(x) 



S(x) = s(x) + W(x) 



(15) 



(16) 



s(x) the classical action and W(x) the quantum action correction function. The structure of (fT5|) enables us to make 
use of the Bakcr-Hausdorff identity, which allows us a formal expansion: 



Q(x) 



exp 



-W{x) 



Q(x) exp 



-W{x) 



Qo(x)--[W(x),Q(x) 



W(x), W(x),Q(x) 



(17) 



= Qq(x) + Q 1 (x)+Q 2 (x) + ... . 
The zeroth order correction function Qq(x) comes from 

Q (x) = Q(x)-l = -e X p[-^2s(x)]J x p'(y)eiq ? [^2s(y)]dy . (18) 
This corresponds to the zeroth iteration of equation (|15p . where the action correction W{x) is chosen zero. 



B. Convergence Properties 

The convergence behaviour of our perturbation series depends crucially on its key quantity, the function Qq(x). 
The whole perturbation series is built on this zeroth order quantum correction function, and thus the behaviour of 
this function will influence the behaviour of the perturbation series. The question we have to ask must be, if there are 
singularities around classical turning points, which plague the WKB-theory. Anticipating the answer, we can state: 
the function Qo(x) is regular in the vincity of classical turning points Xt, if, and only if the potential function U(x) 



4 



is regular there. The examination now follows. 

We want to examine, if the integral, which determines the function Qo(x) posseses upper or lower bounds. We have 

Qo(x) = -exp[i2s(x)] J p'(y)exp[i2s(y)]dy . (19) 

Since we are interested in the behaviour of Qo(x) in the vincity of a classical turning point, we expand the integrand 
around the classical turning point: 



p'(x) = - 



U'{x) 



U'(x t ) 



2^/k* - U(x) ' 2y/-U'(x t )(x-x t ) 



s(x) = I ^ 2 - U(x)dx sa -y/-U'(x t )(x - xtf' 2 



(20) 



This approximation is suitable, because the function p'(x), which is the origin of the turning point singularity is 
approximately zero almost everywhere, except the turning point region. Thus, the main contribution to the integral 
comes directely from the turning point region. Inserting the approximate quantities and performing the integral, we 
get the following estimate 



31/3 



"4~ 




U'{x t ) 


3 




UV*{x t ) 



(21) 



We clearly see, that near the turning point, the behaviour of Qo(x) depends directly on the behaviour of the potential 
function U(x). When U(x) is regular at xt, which is the case for all potentials, that locally can be represented by 
powers x a , a > 0, the quantum correction function itself is regular and finite. When U(x) is unbounded near x t , 
which is the case for potentials, that locally can be represented by inverse powers x~ a , a > 0, then the quantum 
correction function will diverge, too. 

From a more mahematical point of view, Mandelzveig [Tij has shown, that expressions like (|15p converge quadrat- 
ically. Physically spoken, quadratical convergence means, that the zeroth iterate Qo(x) is already very close to the 
exact result Q(x), such that qualitatively the whole quantum mechanical behaviour of a system is to be decribed by 
the zeroth iterate Qo(x). 



C. The Order of the Pertubing Terms 

From the perturbation series (fTT|) , the order in ft is not easily visible. We have seen in (jTTJ) , that the full quantum 
momentum can be expressed as a series 

00 

P(x)=p(x) + ^Q li (x) ■ ( 22 ) 

On the other hand, it is known, that the full quantum momentum can be expanded in power of ft, when ft is small 
compared to the classical action function: 

00 

P(x) = p{x) + ^tfP v {x) . (23) 
i/=i 

The crucial question is now, how these two series are related to each other. To get more insight into this question, 
we start with the selfconsistent expression (jT5j) for Q(x). Due to selfconsistency, the correction function Q(x) is a 
function of all orders of ft. To see this more clearly, we formally can put the indices fi, v together and obtain 

00 00 00 

fi—O fi—0 V— 1 

00 

QM = J2 nu Q^( x ) ■ (24) 

This means, that every order of quantum correction functions carries all orders of ft, as stated above. This astounding 
fact will be demonstrated in the next subsection, where the power series in ft is calculated directely from the quantum 
correction functions. 
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D. Derivation of the WKB-Regime 



After these formal preparations, we will calculate the power series in h. One of the most interesting results of this 
calculation will be a reason for the bad convergence of the WKB-series. 

The WKB-theory follows from the assumption, that the solutions of the Hamilton- Jacobi equation can be represented 
as power series in h: 



S(x) 
P(x) 



s(x) 
p{x) 



J2h u S v (x) 

oo 



(25) 



u=l 



When these ansatze are inserted into the Hamilton- Jacobi equation (0[5]), there arises an infinite hierarchy of equa- 
tions, which can only be solved by successive approximation. Breaking this infinite hierarchy after the first order 
yields the familiar WKB result for the wave-function. But the expansion in powers of h is plagued by singularities 
around classical turning points, and the strength of these singularities increases by every order in h. This is not the 
case, when the series expansion is derived from the quantum correction functions, because their expansion in powers 
of h will create residual terms of higher order in H, which will catch the divergencies and result finite expressions as 
a whole. 

Technically, the expansion can be set up as follows. By use of the auxiliary identity 1 = p(x)/p(x), we get for the 
zeroth order 



Qo(x) 



exp 
exp 



--2s(x) 
-i2s(x) 



p'(y) exp 
'p(v) 



- n 2s(y) 



p'(y) 



p(y) 



exp 



dy 



- h 2 S (y) 



dy 



Performing an integration by parts, one gets 

_ . . i p'(x) i 
Qo(x) = h 2rtx) +h 2 eXP 



-~2s{x) 



dz\ P - 
P 



(z) exp 



-Mz) 



Using the auxiliary identity once again and integrating by parts yields 



Qo{x) 



. i p'{x) 
l 2~pjx) 



h2 l_p(x)p"{x) - {p'{x)f 



p 3 (x) 



+ Ro,a( x ) 



(26) 



(27) 



(28) 



We observe the following facts: the term proportional to H is precisely the first order series expansion following from 
the WKB-theory; the following term is of order ti 2 , and does not contribute to the first order in h. The residual 
term, denoted by i? M= o,j/=3(x) carries only powers of h equal or larger v = 3. Such a residual term is not known in 
WKB-theory, it is lost by initially assuming the power expansion in terms of h. 

Up to now, we have shown, that the zeroth order quantum correction leads to the correct WKB-term first order in h, 



HPtix) =Q ,i(x) 



. i p'(x) 



The second order expansion in WKB-theory gives 



h z P 2 (x) = -hf- 



lp(x)p"(x)-(p'(x)f 



p 3 (x) 



1 {p'{x)f 



p 3 (x) 



(29) 



(30) 



The first term has already been obtained by the second integration by parts of Qq(x). Further integrations by parts 
of -Ro.3(a0 in the same manner would not lead to the required second term, it can only be calculated by expanding 
Qi(x) into powers of h. The perturbational expression for Qi(x) as calculated above (|17p is 



Qi(*) = -fc[/ Q(z)dz,Q(x)} = -^[h Q Q ,i(z)dz,CXx)] 



(31) 
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The first integration by parts cancels, such that only term of order 0(h 2 ) remain, 

Qi(x) = -h l -e-^ 2s ^ J\z^Q ,i(z)e iMz) -\f ' Qo,i{z)dz,Ro, 3 (x)] 

= -ft 2 7^7rGo,i(aO + RM -[C Q Q ,i{z)dz,n^{x)] . (32) 

4p z (x) J 
The first order in ft is already known from above, Qo,i(x) = | p p ^j , leading to the first term of the desired WKB 



expression: 



With that, we obtain 



Qi(x) = Q 1 , 2 (x) + ... = ft 2 iM^ + ... . (33 ) 

8 jr(x) 



P\{x) = Qo.i(x) , 
P*(x) = Q , 2 (x) + Qi,2(x) 



(34) 



The residual terms R^ v {x) und R^x) carry only higher orders in ft. 

With this, rather cumbersome, procedure we could show, how all orders of the WKB-series can be obtained, at least 
in principle. It further clarifies the fidelity of the quantum correction functions: Qo{x) is exact to the first order in ft, 
Qo(x) + Qi(x) is exact to second order in ft and so on. Keeping the WKB-form, we can write the quantum mechanical 
momentum as 

N oc 00 

p{x)= P {x)+Y J wPv(x) + Y. E n ^ x ) > ( 35 ) 

v=l fj,=0v=N+l 

where in 1Z(x) we have subsummized all residual terms. The number N gives the order in ft. The residual terms can not 
be obtained when the WKB-series expansion is set initially. The residual terms result from partially integrating the 
nonsingular expression for the quantum correction function Q(x), but a nonsingular expression will remain nonsingular. 
The procedure of partial integration splits the nonsingular quantum correction function into a singular WKB-term 
and a singular residual term. Both singularities kept together will give a nonsingular function, but omitting the 
residual terms, as WKB-theory does, will give a badly converging, singular power series in ft. 

III. THE CONSTRUCTION OF THE WAVE-FUNCTION 

In the introduction section we have seen, that in our theory the wave-function of a system can be written in terms 
of superponed exponentials, or, if one prefers, in terms of trigonometric functions. Depending on the problem, both 
representations have advantages and disadvantages. The key to the construction of the wave-function is the action 
integral, that depends on the direction of integration. In a two turning-point system, we have the choice to integrate 
from the left (xi) to the right (x 2 ) or vice versa. We have reserved a whole section to this problem, because its 
solution is not as easy as it might look at first. 

We begin with the wave-function, that is obtained by integrating from the left to the right, x\ — ► x 2 . We find, in 
trigonometric fashion, 



^(x) = cxp 



r-X "I r r-X 1 

- QQ{x')dx' cos / (p(x') + $tQ(x'))dx' --^ 



(36) 



The arrow-index of the reflection phase indicates, that the direction of integration here was chosen as x\ — > x 2 . To 
construct this wave-function, we have made the following choice for the quantum momentum function 

P{x) = p(x) + 3?Q(x) + iQQ(x) . (37) 

This choice is of somewhat arbitrary character, we could also have choosen the complex conjugate of P(x) or —P(x) 
or any other possible combination of signs and complex conjugates. But the choice, once made, gives the system the 
property of a direction of integration, all other directions of integration have to be consistent with. This becomes 
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clear, when we now construct the wave- function for the direction of integration x\ <— x 2 . The easiest way to do this 
is to introduce a parity transformation by the operator Tp, that performs an axialsymmetric transformation. The 
position of this axis will be the minimum of the potential - the potential itself may be of unsymmetric shape, the 
minimum is the only condition - that any two turning-point system will possess. We want to know what a wave, that 
starts at x\ would look like, when it starts at x 2 : 



T P exp 



i {p(x') + nQ(x')+i%Q(x'))dx' 



exp 



i { P (x')+$lQ{x)-i%Q(x'))dx' 



(38) 



The sign of the imaginary momentum function has changed. This is due to the fact, that the real parts of the quantum 
momentum are even, the imaginary part is an odd function with respect to the minimum of the potential. This can 
be deduced from (|15[) by inspection. Thus, when a parity transformation with respect to this minimum is applied, 
the sign of the imaginary part must be changed. Without this fact, a consistent construction of a wave-function for 
the direction of integration x\ <— X2 is not possible. We finally obtain 



^>(x) = exp 



exp 



^sQ(x')dx' 



%Q(x')dx' 



(p(x') + "RQ(x')) dx' - 



{p(x') + diQ(x')) dx' 



(39) 



IV. SOFT BOUNDARY PROBLEMS AND QUANTIZATION 



In the framework of Hamilton- Jacobi theory, soft boundary problems arise, because the space is separated into 
classically allowed and classically forbidden regions. By solving the appropriate Hamilton- Jacobi equations, a wave- 
function for each of these regions is known. The boundary value problem is to match all these wave-functions, to 
obtain a wave-function ^ defined on the whole space. This boundary value problem is also known as the connection 
problem, see e. g. In the Schrodinger point of view the connection problem 



*i(:c) <-> V 2 (x) <-> # 3 (a:) 



(40) 



is a formidable task and will, except some few solvable examples, remain unsovled in general. In the Hamilton- Jacobi 
point of view, where we can express the wave-functions as exponentials, the connection problem becomes easily 
treatable on the level of trigonometric and exponential functions: 



exp [— Si(x)] <-* exp [— SJS^a;)] cos 



KS 2 (x) - | 



exp [-S 3 (x)] 



(41) 



The reflection phase 4> solves the connection problem and with <fi the wave-function can be defined on the whole space. 

The points where the wave-function of different spatial regions have to be matched are the classical turning points. 
The boundary value problem, e.g. 



d_ 

dx 



In [^ 2 {x) 



= ^M*3(*) 



(42) 



is called a soft boundary problem, because it is neither a Dirichlet nor a Neumann boundary problem, but a mixture 
of both. Writing (|42l) in the spirit of (|4"Tj) , we obtain 



%Q 2 {x 2 ) - $Q 2 (x 2 ) tan 



f * (p 2 (x')+$tQ 2 (x'))dx' -tzL 

J xi ^ 



= -Qa(x 2 ) 



From this, we obtain the reflection phase 

^Q 2 (x 2 ) - Qz{x 2 ) 



2 arctan 



MQ 2 (x 2 ) 



For the reverse direction, we find similary 

"3fQ 2 (xi) + Qi(x 1 ) 



!><_ = 2 arctan 



( P2 (x')+nQ 2 (x'))dx' 



(p 2 (x')+3tQ 2 (x'))dx' 



(43) 



(44) 



(45) 



8 



The reflection phases, by their dependence on the turning points, are functions of energy. 

Knowing the reflection phases, we can easily derive the quantization law by the argument, that both directions of 
integration must give the same wave- function. The phases of this wave-functions may only differ by mr. From this, 
we get 



r%2 

I nP 2 {x')dx' 

J Xl 



(46) 



Note that by (|4"H H5|) the action-integrals cancel each other in the sum cf)^ + <j>^_ . To make contact with the usual 
Bohr-Sommerfeld rule, we rewrite (146|) and introduce the Maslov-index (j,^ 



'p 2 (x')dx' = n(n+^) , 



7T 



■1-2 



H = - — - / MQ2(x')dx' . (47) 



With (|4T|) the Maslov-index of a physical system is directely calculable by knowledge of the quantum correction 
function. This is a very valuable result, because quantization in the Hamilton- Jacobi point of view means just to 
calculate (|4"7T) to obtain the energy-levels. There is no need to orthonormalize the wave- functions to find the energy- 
levels, when calculated by (|4"7| . The theory of Hilbert-spaces ensures, that wave- functions corresponding to eigenvalues 
of the underlying problem are already orthonormal. 

By (1441 145[) . we can examine another important property of the reflection-phases. They are divided into two parts, 
one gauge-dependent, one gauge- invariant, 

^ = f + r, (48) 

Inspection of (|44l |45|) . in combination with (|36l [39)) shows, that the lower bound of integration does not necessarily 
have to be a turning point. Dependent on the direction of integration, which has to be defined at first, the lower 
bound serves as a point of reference and can be chosen arbitrarily as any point lying between the two turning points. 
The choice of the lower bound must leave the wave- function invariant. This is ensured by the gauge-dependent part 
of the reflection-phase, e.g. 

X2 



0^ = 2/ {p 2 {x') + KQ 2 {x'))dx' . (49) 
The gauge-invariant part is thus given by, e.g. 



Zref 



= 2 arctan 



^Q 2 {x 2 ) - Qz{x 2 ) 



VtQ 2 (x 2 ) 



(50) 



V. QUANTUM TRAJECTORIES 



For some reasons, this paper would be the place to set a profound and philosophical discussion on the sense and 
usefulness of quantum trajectories into work, but we think that should be done elsewhere. Instead, we want to 
demonstrate how easily quantum trajectories are available within the framework of the stationary Hamilton- Jacobi 
theory. The topic of quantum trajectories is often discussed by introducing a quantum potential and a quantum force. 
We refuse the introduction of those quantities, because in our opinion they are somewhat misleading. The so called 
quantum potential is not a potential function as classical potentials are, but depends on the real part of the reduced 
quantum mechanical action function, so does the quantum force. From this point of view, one anyway has to solve 
the Hamilton-Jacobi equation at first and calculate the quantum potential afterwards, which makes the introduction 
of a quantum potential gratuitous (for the here discussed stationary theory). There is another crucial point, that 
leads us to reject the concept of a quantum potential, because, as we will show below, the classical relation between 
momentum and velocity does not hold within the framework of quantum mechanical Hamilton-Jacobi theory. The 
trajectory has thus to be derived from the quantum mechanical action function directely. From this trajectory, by 
time derivatives, it would in principal be possible to construct a Lagrangian or Newton's equation of motion, but such 
attempts are, again, gratuitous to our regards. 
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The use of the stationary Hamilton-Jacobi equation describes isoenergetic dynamics, meaning that motion takes 
place on a energy-shell in phase-space. The presence of quantum mechanical behaviour has no effect on the energy-shell 
itself. This can be demonstrated by writing down the the classical and quantum mechanical Hamiltonians 

2 

n cl = %- + u(x) , 

2m 

P 2 ih d 

H qm = JLp + u(x) , (51) 



and comparing them 



2m 2m dx 



W qm -W c i = 5W . (52) 



If in (|52[) STL is zero, the energy of quantum and classical motion is constrained to the same energy-shell, quantum 
and classical motion are isoenergetic to each other. The employment of quasilincarization P = p + Q naturally fullfills 
STL = because 

Hqm = TL C \ , (53) 

= 2 P Q + Q 2 -ih^- ip+Q) . (54) 

The presence of the quantum momentum corrections thus have no effect on the energy of the system, they only change 
the dynamical flow on the energy-shell. For the dynamics that means, that the time the motion between two points 
takes is changed. In the light of (|S"4")) it becomes clear how misleading the introduction of concepts like quantum 
potentials are, because a genuine potential function should of course change the total energy of the system. 
The presence of a quantum momentum correction leads to a different trajectory, 

Xq(t) = Xd(t) + x qm (t) , (55) 

where the total motion is divised into a purely classical and a purely quantum mechanical part. Naively one would 
think, that from (|55|) it immedeately follows that 



in 



dx q (t) 
dt 



= P (x) + %tQ(x) , (56) 



which, as we will show now, is not true. In classical mechanics, the trajectory equation (|60j) and the velocity - 
momentum relation are equal. We have 

dt = dx ^^^ = (^y 1 = (™y 1 = m . (57 ) 

dE dt \ dE ) \p{x) J m y ' 

In quantum mechanics, the quantum momentum correction breaks this equality, leading to 



It can easily be seen, that the naive assumption (|56|) leads to an obviously wrong result. See [16|, [17| for statements, 
that dig in the same hole. 



The necessary theoretical tool to calculate a trajectory is already known from classical mechanics. As there, we 
start with the action function 

SJS= / dx(p(x) + 9lQ(x))-E(t-to) . (59) 

J a; ro f 

Employing the Jacobi-connection, which is the same as demanding the phase of the wave-function to be stationary, 
it follows 



(60) 
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which we will call the trajectory equation. The use of the Jacobi-connection is a natural consequence of the Hamilton- 
Jacobi theory. Theoretically, it is very well founded and will give the right trajectory for any given action-function. 
Thus, we dwell on a very safe terrain, whereas (|56j) is not more than an assumption. 

It is common knowledge, that a standing wave, e.g. a bound state, has a current density of zero, and thus contains 
no observable motion. As already stated above, this is true in cases where only the wave-function ^ is known. The 
use of the quantum action function S instead allows a decomposition into running waves, which certainly posses a 
dynamical interpretation. This dynamical interpretation allows us to connect the motion of a quantum particle with 
a causal picture. To do so, we have to work out the differences between quantum and classical motion at first. From 
(|60p we get the running time differences between classical and quantum motion moving from an arbitrarily chosen 
point to a turning point and back: 



fX2 

At^ = I dx 



dp(x) d$tP(x) 



v dE dE 



X TQ f 



\ dE dE 



From (|6ip . we obtain the total time-shift between quantum and classical motion 

At = At^ - At«_ = -2 f 2 dx m< ?^ x) . (62) 
Jx re( dE 

The total time-shift is a positive function of energy and a consequence of the isoenergeticity. Its meaning is the 
following: the quantum particle moves with a larger local momentum as the classical particle does and reaches the 
turning point earlier than the classical particle. But for the backward motion the quantum particle is delayed with 
respect to the classical particle, resulting a positive, total time-shift. The delay of the quantum particle at the turning 
point stems from the finite quantum momentum there and reflects the finite probability for entering the potential 
barrier. We emphasize, that a negative sign of the total time-shift (|62[) would result in acausal behaviour, because 
then the quantum particle would leave before it would have arrived, which is not possible. The result (|62[) allows 
us to derive a remarkable property of quantum trajectories, that is, they leave the phase-space integral invariant. 
Comparing the classical and the quantum mechanical phase-space integral, this can easily be seen. If we define T as 
the classical period and 2(Ti + T2) = T, we have 

<j) P (x)dx-ET = I (p(x)+RQ(x))dx-2E{(T 1 + At^) + (T 2 -At^)}^ 

EAt = -~ I $lQ(x)dx . (63) 

From the last line of equation (|63|) equation ((62)) readily follows. Thus, if isoenergeticity is fulfilled, classical and 
quantum mechanical motion are described by the same action hypersurface. With the time-shift (|62p . we now can 
construct the trajectories 

dt = / U dx (x rof ->x 2 ) , 

Jx tat <3E 

dt' = - P K^p-dx' (x ref ^x 2 ) . (64) 



ti+At 



By this procedure one obtains quantum trajectories in the classically allowed region, forming a sum of disconnected 
branches 

00 

x q {t) = Y J (-W l x(t)0{t n+1 -t n - At-\t-t n+1 +t n + At\) . (65) 



The trajectories form disconnected branches, because their shape in the classically forbidden region is not available by 
the trajectory equation (|60|) . but below, when the quantum trajectories of the harmonic oscillator will be discussed, 
we will see how this problem can be solved in a physically reasonable way. 
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VI. APPLICATION 

To demonstrate how quantum mechanics is done from the Hamilton- Jacobi point of view, we will discuss two 
important systems, that are widely used as models for various physical systems: the linear potential and the harmonic 
oscillator potential. Especially, we will work out the differences between the Schrodinger point of view and the 
Hamilton- Jacobi point of view. By the ladder, we will obtain an insight into the quantum mechanical behaviour, 
which can not be gained by the wave-function alone. 

A. The Linear Potential 

We consider a linear potential U(x) = —Fx. By this example, we will focus on the improvement, that is gained by 
the use of the zeroth order quantum correction function Qo(x) in contrast to the usual WKB-approximation. From 
the potential, we obtain local the classical momentum and the local classical action as 



p(x) = ^k 2 + fx , 
s(x) = §j(k 2 + fx) 
2mE _ 



3/2 

2mF 



h 2 ' J h 2 

The zeroth order quantum corrections function in this case can be calculated analytically giving 



(66) 
(67) 

(68) 



h{x) =(-i)-* 



6 



exp[— i2s(x)]Ti/ 3 [— i2s(x)] 



(69) 



During the following, we set / = 1 and k = 0, because the system is translationally invariant on the energy-axis. 
In figure l|T|). we have plotted the phase-space functions of the linear potential. The quantum correction function 




FIG. 1: Classically allowed phase-space of the linear potential, showing the corrected momentum 3?P(a) = p{x) + KQo( x ') (full line), the 
quantum correction function 3?Qo(a;) (dashed line) and the classical momentum p(x) (dotted line). The effect of the zeroth order quantum 
correction can clearly be distinguished. 



3lQo(x), far away from the classical turning point, here x% = 0, is negligible and the classical momentum p(x) is in 
charge. Approaching the turning point region, the quantum correction function ^RQq(x) starts to grow, describing the 
importance of quantum mechanical behaviour in the vincity of the turning point. At the turning point itself, the real- 
part of the corrected momentum 5RP(cc) is finite, while the classical momentum p(x) is zero as usual. This allows only 
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one interpretation: for a classical particle, the potential represents an inpenetrable barrier, and reaching the classical 
turning point leads undoubtedly to a change of the direction of motion. For a semiclassical particle, described by 
$lP(x) — p(x) + $lQo(x), the classical picture is not true anymore. The corrected momentum 9?P(x) is finite at the 
classical turning point, meaning that the potential barrier becomes penetrable with some probability. For a classical 
particle, the potential barrier represents a hard wall, but for a semiclassical particle, the potential barrier becomes 
soft. The softness of the potential barrier in the semiclassical case is decribed by the finite momentum functions 
and by the finite probability to find a particle inside the barrier. We want to emphasize here, that in the WKB- 
approach the potential barrier remains hard, because of the turning point singularity. The proper connection between 
classically allowed and classically forbidden region in the WKB-case can only be established by the knowledge of the 
exact solution of the Schrodinger equation, leading to the reflection-phase cj> = ir/2. By the use of quantum correction 
functions, the turning point singularity is overcome. The linear potential serves as a model for the one turning point 




FIG. 2: Comparison of the wave- function obtained by the zeroth order quantum correction function (full line) and the usual WKB- 
approximation (dashed line). Already the zeroth order quantum corrections removes the WKB-turning point singularity completely and 
is, up to the amplitude, identical with the exact solution of the corresponding Schrodinger equation, the Airy function. 



case. Taking the turning point X\ as the turning point of the system, the wave-function can be constructed by 



&(x) = exp 



(x')dx' 



cos 


-/ 




J X 



(70) 



The second turning point xi must be set to infinity. This becomes clear when we remember that (|39[) was derived 
from the symmetry properties of a two turning point system: the particle comes in from infinity until x\ is reached, 
where it is reflected back to infinity - and at infinity it is reflected again into direction of X\. Thus, a one turning point 
system follows from a two turning point system by setting one of the two turning points to infinity. This procedure is 
necessary, because it is the only way to take into account the symmetry properties of the quantum correction function 
as discussed above. Together with (1451) and the symmetry property of *<sQo(x) we obtain the wave-function 



^(x) = exp 



QQ (x')dx' 



cos 



(p( x ')+^Q ( x '))dx'-^f- 



(71) 



Expression (|7Tj) perfectely demonstrates how the gauge-dependent part of the reflection-phase cancels, while the 
gauge-invariant part <^> inv remains present. As can be seen in fig. the invariant phase </> mv = 2ir/3 is not the WKB 
reflection-phase n/2 obtained by the asymptotical expansion of the Airy- function. The asymptotical reflection-phase 
follows from the limit 



1 J Xl 



lim / UQ (x')dx' 

x—>oo 



7T 

4 



(72) 
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FIG. 3: Illustration how the correct asymptotic reflection-phase is approached by the interplay ol the real part of the quantum action 
correction and the kinematical phase. 



The harmonic oscillator is surely the most considered and best understood quantum system. This fact makes 
it a very interesting test object for the viewpoint on quantum mechanics we haves developed here, because by the 
harmonic oscillator, we can demonstrate at best the differences between Schrodinger and Hamilton-Jacobi quantum 
mechanics. 

We will work in suitably scaled quantities 



For simplicity, but without loss of generality, we have set 17 = 1 in all numerical calculations. In figure (f4]) we have 
plotted the quantum mechanical and the classical phase-space of the harmonic oscillator for the ground state energy 
(left picture). Again we see the typical effect of the quantum correction functions, breaking up the closed classical 
phase-space, thus making the potential barrier soft and penetrable. The effect of the quantum correction functions 
is largest for the ground state, as one would naturally expect, because the ground state is close to the anticlassical 
limit of the system. For the tenth excited state, the right picture in ([4|) shows that the influence of Q(x) already has 
declined, the system has already entered its semiclassical limit and quantum corrections become negligible. 

We will not present pictures of oscillator wave- functions here, because we think that up to now the reader might 
be convinced, that the wave- functions determined by (|36p or (139p are, up to the normalization amplitude, identical 
with the oscillator wave-functions determined by the Schrodinger equation. Instead we will focus on the calculation 
of quantities, that are not available within the Schrodinger equation. In figure ([5]) we have calculated the quantiza- 
tion curve and the Maslov-index according the formulae (|46[ I47[) . The quantization curve determines the oscillator 
eigenstates by knowledge of the classical momentum and its quantum corrections. Although the wave-function is 
implicitly known when the reflection phases (|44| I45|) are known, there is no need to employ them for determining 
the spectrum. The knowledge of the quantum correction functions is sufficient for a calculation from first principles. 
The Maslov-index imposes a tiny oscillating structure onto the quantization curve. The nodes of these oscillations lie 
perfectely at 1/2, coincidence with the bound states and recover the Bohr-Sommerfeld rule. This behaviour of the 
Maslov-index was first calculated in [l5| for the half-sided harmonic oscillator and by a different method. 



B. 



The Harmonic Oscillator 




(73) 
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FIG. 4: Comparison of the phase-space of the harmonic oscillator for the full quantum momentum 5RP(x) = p(x) + ^RQ(x) (full line) 
and the classical momentum p(x) for the ground state (left) and the tenth excited state (right). Due to the soft boundary conditions, the 
quantum mechanical phase-space is an open ellipse, while the classical phase-space is the well-known closed ellipse. 



C. Trajectories for the Harmonic Oscillator 



For the harmonic oscillator, for further details see part I of this paper, we now will discuss the topic of quantum 
trajectories in detail. We will calculate and discuss trajectories for the ground state and the tenth excited state to 
illustrate the influence of the quantum corrections functions near the anticlassical and in the semiclassical regime. 

By (|62p the time-shift between quantum and classical motion, sec figure ^ is given as 



The lower bound in (|74p is zero, because we want to calculate the trajectory of a particle starting its journey from 
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FIG. 5: Quantization curve and for the harmonic oscillator in units of ir. 



the minimum of the oscillator potential x — 0. We could have chosen any other point, the choice of the point of 
reference has no influence on the shape of the quantum trajectory. From the trajectory equation (|60p . we get 



x q (t) 



n=0 



-l) n x(t)9(t n 



+i 



t, 



At - \t - t n+1 + t n + At\) 



(75) 



where the index n denotes the multiple branches obtained by shifting the time by At in the classically allowed region 
for the forward and backward motion, respectively. By the trajctory equation (|60|) only the quantum motion between 
the two turning points can be determined. The cutoff at the turning points is the cause the multiple branches. If we 
compare (T75|) with the classical trajectory 



x c \(t) 



shx[Clt] 



(76) 
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FIG. 6: Total time-shift At as function of ft. For high energies, the total time-shift tends to zero, because the system approaches its 
classical limit, where quantum and classical trajectories do not differ anymore. 

and unse the fact that quantum and classical motion are described by the same frequency, see (|63|) . we are lead to 
make the following ansatz for a quantum trajectory 



The limit in (177j) requires, that the quantum amplitude approaches its classical value in the classical limit. The fact 
that all quantum effects are now located in an amplitude-function reminds us of the path-integral, where the same 
phenomenon occurs: phase-quantities are given by the classical action function, but the time-evolution amplitude is 
quantum mechanical. However, by comparison between (|751 177[) . we can establish a continous quantum trajectory 
without cutoffs at classical turning points by analytical continuation beyond the classical turning points, illustrated 
in figure ([7]). The reader might object this procedure by the argument, that by analytical continuation in this case 
a kind of information is gained, that does not exist. We agree with this argument to some extent and will therefor 
call the trajectories in the classically forbidden region apparent trajectories. The present situation is similar to the 
quantum time-shifts in [18j |. where it was found that apparent trajectories enter the classically forbidden region. But 
the difference between there to present is, that in [l8j the potential was a singular one - and no particle, classical or 
quantum can cross a singularity. But the harmonic oscillator is regular and by the quantum correction functions, the 
turning point region has become regular, too. Further know we from above, that in the quantum case the potential 
becomes soft and therefor penetrable. All this put together means the following: for the harmonic oscillator, the 
reflection probability at a turning point is \R(k n )\ 2 — 1 for all energies, thus the particle comes back in any case. The 
particle reaches the classical turning point with finite momentum, which imposes a time-shift in comparison to the 
classical particle and a phase-shift on the wave-function as a result of the soft boundary condition to be fullfilled. The 
phase-shift of the wave-function and the time-shift of the motion are two manifestation of the same phenomenon, the 
relative softness of the potential barrier in the quantum mechanical limit. All these facts allow only one interpretation, 
the particle as it has arrived at the turning point with finite momentum, goes beyond the turning point and enters the 
classically forbidden region, appears to be decelerated by the soft barrier, comes to rest and finally moves backward to 
be expelled again into the classically allowed region. The only explanation thus is the softness-effect of the potential 
barrier. We emphasize, that this softness has nothing to do with sticking or absorptive reactions, it is a pure quantum 
mechanical effect. For energies approaching the classical limit, the barrier hardens as is described by the limiting 
process in ((77)) . 




(77) 
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FIG. 7: Comparison of quantum trajectories (full and dotted lines) and classical trajectories (dashed lines). The full lines are the branches 
obtained by the trajectory equation H60I I. The total time-shift At is the amount of time the quantum particle spends in the turning point 
region. It can be clearly seen, how the limiting process from l|77| l works. 



VII. SUMMARY AND CONCLUSION 

We have shown how quantum mechanical problems can be solved within the framework of the Hamilton- Jacobi 
equation of quantum mechanics, where the central quantity of interest is the quantum mechanical action function 
instead of the wave-function. To our regards, Hamilton- Jacobi theory and Schrodinger theory are fully equivalent to 
each other, they describe the same physics and differ only in their descriptive manner. 

In contrast to the well-known Madelung-de-Broglie-Bohm theory, we do not separate real and imaginary part of 
the action function initially. A very valuable feature of Hamilton-Jacobi theory is, that by use of the quantum 
action function, the wave-function can uniquely be decomposed into running waves in every point in space, regardless 
whether the wave-functions describes bound or scattering states. In the Schrodinger point of view, this is possible 
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only asymptotically. 

The use of quantum correction functions enopens a deeper insight into the quantum mechanical behaviour of 
systems, especially in the vincity of turning points, where the local quantum mechanical momentum is finite. This 
property overcomes the WKB-singulariy at classical turning points and furthermore clearihes the special behaviour 
of a potential barrier in the quantum mechanical case: the barrier becomes soft and therefor allows the quantum 
particle to enter it. This softness-property of the potential barrier in the quantum mechanical case is modeled by 
the soft boundary conditions the wave-function has to obey - and that can be solved exactly by help of the quantum 
action function - and also explain the finite probability of finding a particle there. 

The soft boundary-value problem is solved by a reflection phase. The reflection phase is the key-quantity to obtain 
the spectrum of a quantum system with help of a quantization curve and allows a direct calculation of Maslov- 
indices. To obtain the spectrum of a quantum system by this method, only action and momentum functions have 
to be known, orthonormalized wave-functions are not necessary. Especially for potentials where exact wave-functions 
are not analytically known, it may be much more simple to calculate the quantization curve instead of searching 
orthonormalized eigenstates with the Schrodinger equation. 

Besides the Hamilton- Jacobi theory allows us to calculate quantities that are not available within the Schrodinger 
point of view, the calculation of quantum trajectories is possible. The quantum trajectories obtained by the trajectory 
equation (|60| are only defined in the classically allowed region, but by analytical continuation, they can be established 
also for the classically forbidden region in a physically reasonable manner. The behaviour of the trajectories again 
mirror the softness-property of the potential barrier in the quantum mechanical case and the finite probability to 
enter the classically forbidden region is all described by a total time-shift At (|^|) . In contrast to de-Broglie-Bohm 
theories, no quantum potential is introduced, and therefor no equation of motion. Moreover, the equation of motion 
in de-Broglie-Bohm theories [!, 0, [I[ and our trajectory equation ([SO)) contradict each other. This result was also 
found by [7, 8j, but was interpreted differentely. Although the trajectories obtained from the quantum action integral 
by use of the Hamilton-Jacobi formalism show a somewhat quirky behaviour, they are consistent with our results 
about the behaviour of the quantum correction function in the vincity of classical turning points and the behaviour of 
the wave-function. Both, quantum correction functions and trajectories illustrate the softness of a potential function 
near the anticlassical limit. Moreover, Hamilton-Jacobi formalism is theoretically very well founded and thus we 
expect ourselves so travel on safe terrain, whereas, as stated in [l6l.ll7| the equivalence between velocity and quantum 
mechanical momentum in de-Broglie-Bohm theories is up to now an unproved assumption. 



APPENDIX A: COMMENT ON QUANTUM POTENTIAL THEORIES 



The work of Floyd [?J interprets stationary Hamilton-Jacobi theory in a very different way as we have done here. 
Although starting from the same key equation @, Floyd developed a completely different theory by separating real 
and imaginary parts of the quantum action function. By the help of a constrained condition similar to (|54p he arrived 
at a nonlinear, second-order equation for the real part of the quantum momentum function 



E = P 2 {x) + U{x) + Q(x) 



Q(x) 



2 



P"(x) 
~P(x) 



P(x) 



(Al) 



where Q(x) denotes a quantum potential. If P(x) is the classical momentum p(x) 7 we obtain from (|A1[) the well-known 
quantality function, that describes the break-down of WKB-theory in the vincity of classical turning points, see [lOj. 
But if we are interested in in the quantum momentum, we are not allowed to simply write 



P{x) = VE-U(x)- Q{x) 



(A2) 



because Q(x) depends on P(x). Moreover, if P{x) has been obtained by (|Alj) . equation (|A2j) reduces to an identity 
P{x) = P{x) which is of no use either. Because of this, we have argued above, that the concept of a quantum potential 
becomes misleading in stationary Hamilton-Jacobi theory. 

Our second objection is, that by transforming ([2]) into a real but second order equation, there occurs a loss of 
information, because for (|Alj) neither a closed solution formula, nor any initial conditions are known. The apparent 
arbitrariness for the choice of initial conditions in (|A1|) lead Floyd to interpret, that there are infinitely many quantum 
microsates, associated with infinitely many possible initial conditions, that map onto the same energy-eigenstate. For 
nonlinear, higher order differential equations, its the crux of matter that no unique solution may exist, instead of 
a large number of different solutions, depending sensitively on initial conditions. But from such a multitude of 
solutions, the physically relevant one must be extracted. In our first order approach, this is naturally done by (jlll 
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[T5| . Moreover, the action deduced from (|A1[) . the real part of our action function, is assumed to be real on the whole 
space. Comparing this with our findings above we can state that this is impossible, because a necessary condition 
for the wave-function to fall off in the classically forbidden region is a purely imaginary action function there. If the 
concept of the quantum potential would be dropped, this could be deduced from (|A1|) . together with the equations 
that relate real and imaginary parts of the action function. 

We close this comment and emphasize again, that to our regards, the Hamilton- Jacobi theory of quantum mechanics 
does not offer a kind of hidden information that is lost in the Schrodinger point of view. The information provided 
by both theories is the same, only their descriptive manner is different. 
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